Theory of Elastic Interaction of the Colloidal Particles in the Nematic Liquid Crystal 

Near One Wall and in the Nematic Cell 
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We apply the method developed in Ref. [S.B. Chernyshuk and B.I.Lev, Phys.Rev.E, 81, 041701 
(2010)] for theoretical investigation of colloidal elastic interactions between axially symmetric parti- 
cles in the confined nematic liquid crystal (NLC) near one wall and in the nematic cell with thickness 
L. Both cases of homeotropic and planar director orientations are considered. Particularly dipole- 
dipole, dipole-quadrupole and quadrupole-quadrupole interactions of the one particle with the wall 
and within the nematic cell are found as well as corresponding two particle elastic interactions. A 
set of new results has been predicted: the effective power of repulsion between two dipole particles 
at height h near the homeotropic wall is reduced gradually from inverse 3 to 5 with an increase of 
dimensionless distance r/h; near the planar wall - the effect of dipole-dipole isotropic attraction is 
predicted for large distances r > rdd = 4.76fa; maps of attraction and repulsion zones are crucially 
changed for all interactions near the planar wall and in the planar cell; one dipole particle in the 
homeotropic nematic cell was found to be shifted by the distance 8 eq from the center of the cell 
independent of the thickness L of the cell. The proposed theory fits very well with experimental 
data for the confinement effect of elastic interaction between spheres in the homeotropic cell taken 
from [M.Vilfan et al. Phys.Rev.Lett. 101, 237801, (2008)] in the range 1 -f- 1000/fcT. 



PACS numbers: 61.30.-v.42.70.Df.85.05.Gh, 47.57.J- 



I. INTRODUCTION 

Colloidal particles in nematic liquid crystals (NLC) 
have attracted a large amount of research interest over 
the last few years. Anisotropic properties of the host 
fluid - liquid crystal - give rise to a new class of colloidal 
anisotropic interactions that never occur in isotropic 
hosts. The world of anisotropic liquid crystal colloids 
is much more varied than that of isotropic liquids. 

Experimental study of anisotropic colloidal interac- 
tions in the bulk NLC has been made in These 
interactions result in different structures of colloidal par- 
ticles, such as linear chains along the director field, for 
particles with dipole symmetry of the director, [l|, Q 
and inclined chains with respect to the director for 
quadrupole particles [2]- [5]. Colloidal particles sus- 
pended at the nematic-air interface form 2D hexagonal 
structures 0, [13]. Quasi two-dimensional nematic col- 
loids in thin nematic cells form a rich variety of 2D crys- 
tals by using laser tweezers. There are 2D hexagonal 
quadrupole crystals [HI, [12], anti- ferroelectric dipole type 
2D crystals [ill, EH and mixed 2D crystals [3] sand- 
wiched between cell walls. Levitation effect of dipole par- 
ticles in planar cells was studied in . Long ranged elas- 
tic interactions between colloids have been experimen- 
tally found to be exponentially screened (confinement ef- 
fect) in the nematic cell across distances compared to the 
cell thickness L [l6[ . Experimental results are reproduced 
by using the Landau - de Gennes free-energy numerical 
minimization approach [Il|.[l6j-[l8j as well as molecular 
dynamics [26]. 

Theoretical understanding of the matter in the bulk 
NLC is based on the multipole expansion of the direc- 



tor field and has deep electrostatic analogies (2(|- [25[. In 
spite of some differences in these approaches, only one 
of them [20| gives an exact quantitative result which has 
been proven experimentally. Authors of [|| measured di- 
rectly the dipole-dipole interaction of iron spherical par- 
ticles in a magnetic field and found it to be in accor- 
dance with [201 ] . within a few percents accuracy for a 
dipole term. Authors of 0] and Q have measured exper- 
imentally dipole-dipole interaction and found it to be in 
accordance with [2fJ] within about 10% accuracy. This 
allows to justify main assumptions of [20l ] for spherical 
particles in infinite nematic liquid crystal. 

In spite of some understanding of the matter in the 
bulk NLC, there was no theoretical approach which was 
able to make exact quantitative predictions for nematic 
colloids in the confined liquid crystal, though almost all 
experiments in the nematic liquid crystal field are being 
conducted exactly in the confined volumes like walls of 
interface, cells etc. In papers [27, 28], authors explained 
qualitatively screening effects within the coat approach 
[23| developed for the case of the homeotropic cell. But 
the coat approach can not give exact quantitative results 
for potentials as the parameters of the coat remain un- 
known. 

Recently, a proposal was made to take another ap- 
proach [30] for the description of colloidal particles in 
the confined NLC, which may be a generalization of the 
method [2(| for confined NLC. Using this approach, it is 
possible to find all long-range asymptotic behavior of the 
colloids in confined NLC and to make exact quantitative 
predictions which can be tested and compared with the 
experiment. 

In this paper we apply the proposed Green function 



2 



method [3(| for quantitative description of the elastic col- 
loidal interactions between axially symmetrical colloidal 
particles near one wall and in the nematic cell. We pre- 
dict many new effects that fit very well with experimental 
data for the confinement effect of elastic interaction be- 
tween spheres in the homcotropic cell taken from Il6| in 
the range 1 -j- 1000/cT. 

The outline of the paper is the following : In Sec. 
II we talk about an exact approach for colloidal parti- 
cles in nematic liquid crystals and topological defects. 
Sec. Ill presents the general Green function method for 
the description of colloidal particles in confined NLC, 
Sec. IV presents the energy of one particle near one 
wall with both homeotropic and planar anchoring con- 
ditions, in Sec. V and Sec. VI we present the energy 
of elastic interaction between two particles near one wall 
with homeotropic and planar anchoring conditions, re- 
spectively. Sec. VII presents the energy of one particle 
in the NLC cell with a thickness L with homeotropic 
and planar anchoring conditions. In Sec. VIII we find 
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where n is the unit vector pointing average orientation of 
the long axes of LC molecules (we don't take into account 
K24 and if 13 terms for simplicity, though they exist in 
the full deformation energy. Full energy derivation of 
NLC with K 2 4 and Ki 3 terms from microscopic theory 
was obtained in (3lj|). In the one constant approximation 
K = Ku the total bulk deformation energy has the form: 

F bu i k = y J d 3 x [(V • n) 2 + (V x n) 2 ] (2) 

The magnitude of the elastic constant is K s» IQpN . 

On the surface of the particle, LC molecules have the 
tendency to lie either perpendicular or parallel to the 
surface. The orientation depends on the coating of the 
surface. The resulting surface energy can be written in 
the Rapini-Popula form: 

F SU rfa.ce = W j> da{nvf (3) 

with v normal vector to the surface, W < corresponds 
to the normal orientation of LC molecues (homeotropic 
anchoring) and W > corresponds to the parallel orien- 
tation of LC molecues (planar anchoring) . Typical scale 
of the anchoring coefficient is W = 5 • 10~ 5 J/m 2 . Just 
this surface term is the origin of bulk director deforma- 
tions. For small particles with size R < r c = ^ ss 0.2/im 
deformations are small but when the size of the particle 
exceeds R > r c = |£ ~ 0.2/im , surface energy plays a 
dominant role and the director field follows the surface. 



interactions between two particles within a homeotropic 
nematic cell and in Sec. IX interactions between two par- 
ticles within a planar nematic cell are found. In Sec. X 
we make an analysis of the obtained results and discuss 
them in comparison with the results of other authors . 
And in Sec. XI we make our conclusions. 



II. THE EXACT APPROACH FOR COLLOIDAL 
PARTICLES IN NEMATIC LIQUID CRYSTALS. 
TOPOLOGICAL DEFECTS. 

In this section we show the exact general formulation of 
the problem for a system of colloidal particles in nematic 
liquid crystal. 

Let's consider first one colloidal particle of the size 
0.1 -f- 10/UTO inside the unlimited nematic liquid crystal. 
The free energy of the system consists of the deforma- 
tional energy of the LC and surface energy of the par- 
ticle. Dcformational energy can be written in the well 
known Frank form: 



(1) 
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Naturally, topological defects appear near the particle. 
Fig.l demonstrates examples of topological defects near 
the spherical particle [I?} • 




a) b) 



Figure 1. Topological defects, a) Dipole configuration. Point 
defect - hyperbolic hedgehog, b) Saturn ring quadrupole con- 
figuration with c > and c) Two point defects - boojums at 
the poles of the sphere, quadrupole configuration with c < 0. 

Symmetry of the defects can be either dipole (case a.) 
or quadrupole (cases b. and c. ) for large distances. 
Such topological defects can not be treated precisely ana- 
lytically, but, rather, with the help of approximate varia- 
tional ansatzes or with the help of numerical simulations. 

If we have a system of many colloidal particles we need 
to replace with: 

Fsurface = ^ W i f d(j{nuf (4) 

where i enumerate particles and to find the global min- 
imum of the functional F = Ft> u ik + F sur f ace . Moreover 



F bu i k = \j&* {^i( V ' n ) 2 + #a(n • V x n) 2 + K 2 (n x V x n) 2 } 
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if we confine the liquid crystal with some bounding walls 
we need to add surface energy at each interface: 

F waUs = V / W^da^vf (5) 

where E 7 are confining walls. Then we need to find the 
global minimum of the total free energy F — Fj,^ + 
Fsurface + F wa u s . Since the problem cannot be solved 
precisely, even for one particle, without confining walls; 
it is much more difficult to solve it for a system of many 
particles that are found in the presence of bounding walls. 

Therefore it is necessary to introduce some another 
approach which can effectively describe the system of 
colloidal particles in the confined nematic liquid crys- 
tal. This analytical approach we will develop in the next 
section. 



III. GENERAL GREEN FUNCTION 
APPROACH FOR THE DESCRIPTION OF 
COLLOIDAL PARTICLES IN CONFINED NLC 

Here we will describe, briefly, the method proposed 
in [30| . Consider axially symmetric particle of the mi- 
cron and sub-micron size which may carry topological 
defects such as hyperbolic hedgehog, disclination ring or 
boojums. In the absence of the particle ground, the non- 
deformed state of NLC is the orientation of the director 
n||z,n = (0,0, 1). The immersed particle induces defor- 
mations of the director in the perpendicular directions 
= x,y. In close vicinity to the particle, strong 
deformations and topological defects arise, but beyond 
them, deformations become small. In the paper [23| this 
area with strong deformation and defects was called the 



which brings Euler-Lagrange equations: 

An M = An [9 M P(x) - 9*9 M C(x)] (9) 

where -P(x) and C(x) are dipole and quadrupole mo- 
ment densities, // = x, y and repeated \i means sum- 
mation on x and y like <9 M n M = d x n x + d y n y . For 
the infinite space the solution has the known form: 
n„(x) = /d 3 x'^[-^P(xO + a^C7(xO]. If we 
consider P(x) = pS(x) and C(x) = c5(x) this really 
brings n x {v) = pf, + 3cf§,n y (r) = p% + 3c^. This 
means that effective functional ((HJ correctly describes 
the interaction between particle and LC via linear terms 
— 47rP(x)9 M n M — 47rC(x)9 z 9 /J n (tt . So it can be used, as 
well, for the description of colloidal particles in confined 
nematic liquid crustals. 



coat. Beyond the coat, the bulk energy of deformation 
may be approximately written in the form: 

F bu i k = —J d 3 x(Vn„f (6) 

with Euler-Lagrange equations of Laplace type: 

A? V = (7) 

Then director field outside the coat in the infinite LC has 
the form n x {r) = p^ + 3cf#, n„(r) = p^ + 3cff with p 
and c being dipole and quadrupole elastic moments (we 
use another notation for c with respect to the c in [201 ] . 
so that our c = |c ). Actually p and c are unknown 
quantities. They can be found only as asymptotics from 
exact solutions or from variational ansatzes. It was found 
in [2(| that p = ao? , c = —(3a 3 with a being the particle 
radius, and for instance a = 2.04, (3 = 0.72 for hyperbolic 
hedgehog configuration (see FigQ] a ) from ansatz [2(| ■ 
Experiment @ gives a — 2.05, (3 — 0.2 ±0.1, experiment 
gives a = 1.88 ±0.18, (3 = 0.52 ±0.12, experiment 
gives a = 2.21 ± 0.04, f3 = 0.497 ± 0.09. Consider that 
we have found p and c and they are fixed. Then we can 
formulate an effective theory which describes particles 
and interactions between them, remarkably well. In other 
words, constants p and c are bridges between the effective 
theory and exact theory with free energy @ and 

In order to find effective energy of the system: par- 
ticle^) + LC it is necessary to introduce some effec- 
tive functional P e // so that it's Euler-Lagrange equations 
would have the above solutions. In the [20J it was found 
that in the one constant approximation with Frank con- 
stant K, the effective functional has the form: 



(8) 
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In the case of confined nematic LC with the boundary 
conditions n M (s) = on the surfaces E (Dirichlet bound- 
ary conditions) the solution of EL equation ([§]) has the 
form: 

n„(x)= / dVG(x,x')ha;P(x')+^C(x')] (10) 

where G is the Green function A x G(x, x') = — 4-7r<5(x— x') 
for x, x' G V (V is the volume of the bulk NLC) and 
G(x, s) = for any s of the bounding surfaces E. The 
mathematical symmetry property G(x, x') — G(x',x) 
can be proved for the Green functions satisfying the 
Dirichtle boundary conditions by means of Green's the- 
orem f32j]. 

Consider N particles in the confined NLC, so that 

p(x) = J2iPi s ( x - x and = J2i c ^( x - x 0- 



Feff = K 



Jd*X 



(Vn M ) 2 



47rP(x)5 M n A1 - 47rG(x)<9 z <9 M n A1 



4 



Then substitution ((TO) into F eff brings: F eff = U self + 
jjinteraction where jjseif = j^. U- elf , here U- elf is the 
interaction of the i-th particle with the bounding sur- 
faces U* 01 * = U^ff + U s Jq + Uqq . In the general case, 
the interaction of the particle with bounding surfaces, 
(self-energy part), takes the form: 



rjself 

u dd 



-2^ 2 5 M a;i/(x l ,xO| Xl , 



(11) 



a) 



b) 



c) 



U 



self 

dQ — 



U. 



self _ 



<><> 



-2TrKc 2 d z d' z dM( Xl ,^\ 



Figure 2. Cones of repulsion and attraction for bulk inter- 
actions from (|13[) . '-' means attraction, '+' means repulsion, 
a) Dipole-dipole interaction b) Dipole-quadrupole interaction 
with greater particle in the center c) Quadrupole-quadrupole 
interaction. 



where G(x, x') = 



-x'| 



#(x,x') and A x #(x,x') = 



(we excluded the divergent part of self energy from 
i ) 

|x-x'| >■ 



Interaction energy u mteractwn = ^ 
is the interaction energy between i 
u;f = U dd + U dQ + U QQ : 



jrint 

and 



Udd = -47rA>p'a Ai a'G(x. i ,x' 



Here Wf 
particles: 



(12) 



U 



dQ 



Uqq 

AttK 



-(9-90cos 2 6» + 105cos 4 6») 



where r is the distance between particles, 6 is the angle 
between r and z. 

Below we will apply the expressions , (fT2]) for par- 
ticular cases of the NLC confined with one wall as well as 
for NLC confined with two parallel walls (nematic cell) 
with homeotropic and planar boundary conditions . We 
shall consider below p — era 2 , c = —(3a 3 with a being the 
particle radius. 



U QQ = -AnKcc'd^d^G^,^) 

Here unprimed quantities are used for particle i and 
primed for particle j. U d d,U d Q,UQQ means dipole- 
dipole, dipole-quadrupole and quadrupole-quadrupole in- 
teractions, respectively. 

Formulas (JTT) and (|T2j represent general expressions 
for the self energy of one particle, (energy of interaction 
with the walls), and interparticle elastic interactions in 
the arbitrary confined NLC with strong anchoring con- 
ditions n M (s) = on the bounding surfaces. 

When bounding surfaces are located far away from the 
particles and their influence can be neglected, we have 
the approximation of unlimited nematic liquid crystal. 
In unlimited nematic i/(x, x') = and G — Go(x, x') = 
. Then self energy of the one particle is zero as 
we do not take into account divergent part of self energy 
from Gn = t — — tt. As well we do not take into account 

u |x— x'| 

actual finite self energy (inside coat region) of the one 
particle from the exact free energy (|2) and (|64[). Elastic 
interaction between particles in the infinite LC then has 
the well known form 12011 : 



Udd 
4ttK 



3cos 



IV. INTERACTION OF THE ONE PARTICLE 
WITH THE WALL 

A. Interaction of the one particle with a 
homeotropic wall 

We choose coordinate system for the wall with 
homeotropic conditions ^Hn^, z — at the wall and 
z > above, so that particles and NLC under the wall 
have z < and heights h = —z (sec FigJ3]). Then Green 
function in this case has the form [32[ : 



(-rwallr l\ _ 
^hornX*-! x ) — 



1 



1 



|a(x) 



(14) 



(x,y, —z) for x = (x, y, z). In other 
is the 'potential' of the mir- 



Here a(x) = 
words H (x, x') — , 

v 5 ' |a(x)— x'| 

ror image of the particle located in the point x simi- 
lar as in electrostatics. Then interaction of the par- 
ticle with the wall consists of three parts [/fr 01 — 



u dd,self 

(HJJ: 



U" 



,self 



+ u, 



QQ,self 



self 

and may be found via 



U. 



iom,wali 
dd,self 



nKa 2 a 4 
2h 3 



(15) 



UdQ 

4nK 



(pc' 



cp) — j-(15cos z 6 



9) (13) 



U. 



hom,wai 
dQ,self 



STrKafla 5 



2h 4 



(16) 
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a) 

I I I I I I I I I I 



I P>0 



A 



b) 

I I I I I I I I I I 



p<0 



to the cell plane (x £ (0, oo)): 



s~lWQ.ll / /S _ 

'-'plan V-"-! *■ J ~~ 



1 



1 



|x-x'| |b(x)-x'| ^ 
Here b(x) = (— x, y, z) for x = (x,y,z). In other words 
H(x,x') = — | b ( x j_ x /| is the 'potential' of the mirror im- 
age of the particle located in the point x. Taking deriva- 
tives of the H (x, x') we come to the interaction energy 
of the particle with the planar wall via (fTTj) : 



2„4 



u. 



dd,self 



3irKa 2 a 



(19) 



Figure 3. Repulsion of the one particle from the wall with 
homeotropic director orientation, a) p-configuration b) h- 
configuration. 



u, 



hom,wall 
QQ,self 



3irK(3 2 a 6 
2ft 5 " 



(17) 



where h is the distance from the particle to the wall, '-' 
corresponds to p > (p-configuration) and corre- 
sponds to p < (/i-configuration, see FiglSJ) ■ Formulas 
(P~5]) - (fT71) a s well, may be obtained from Lubensky ap- 
proach [201 ] just considering that there is an image par- 
ticle at the height h above the wall with opposite dipole 
moment p' = —p and the same quadrupole moment d = c 
and taking into account that interaction energy should be 
divided by two, as there is no real liquid crystal above 
the wall. For instance, energy of interaction between two 
dipole particles is U dd = i7r % pp ' (1 - 3cos 2 6), then tak- 
ing p' = —p = —aa 2 ,6 = 0, r = 2h and divided by two 
we obtain formula (|15[) . However this method cannot 
be applicable for the wall with planar orientation of the 
director. We shall demonstrate this in the subsection 
below. 



B. Interaction of the one particle with planar wall 

We choose coordinate system for the wall with planar 
conditions with z||hoq, axis x looks down < x < oo, 
x = at the wall so that particles and NLC have x > 
and heights h = x (see FigfJJ. 

In order to find the Green function for this case, 
let's turn coordinate system (CS) of the homeotropic 
cell CS hom (x,y,z) round the y axis on tt/2. Then 
we will have CS plan (x, y, z) with transition matrix A: 
x = Ax,x' = Ax.' so that x = z,y — y, z — —x. Then 
Gh om (x,x') = G hom {Ax., Ax!) = Gpia„(x,x). 

Omitting sign ~ we may write the Green function for 
planar cell in the CS plan with n\\z and x perpendicular 



U, 



jjplan^wall p. 

U dQ,self — U 



plan.wall _ 15TrK/3 2 a 6 
QQ,self ~~ 16ft, 5 



(20) 



(21) 



Figure 4. Repulsion of the one particle from the wall with 
planar director orientation. 

We see that these results differ from (pT5 |) - (fT7|) for 
homeotropic wall. This means that planar orientation 
of the director on the wall violates the direct analogy 
between nematostatics and electrostatics, due to the dif- 
ferent symmetry of the planar wall. Thus, interaction of 
the one particle with the planar wall cannot be treated 
using the previous results of [2(| as in the case of the 
homeotropic cell. 



V. INTERACTION BETWEEN TWO 
PARTICLES NEAR ONE HOMEOTROPIC WALL 



If we have the Green function (|T4|) . we can find inter- 
action between two particles using formulas (|12[) . Then 
taking derivatives, brings all necessary potentials of elas- 
tic interaction between two particles near one wall with 
homeotropic conditions: 



6 



jjwall 
dd.hom 

AirK 



fa 



3cos 



(22) 



jjwall 
dQ,hom 

AttK 



— (pc — cp) — — (l5cos 2 9 — 9) + (pc 



cp')^-(15cos 2 § - 9) 



(23) 



jjwall 
u QQ,hom 

AllK 



-^-(9 - 90cos 2 6> + 105cos 4 6>) + — (9 - 90cos 2 9 + 105cos 4 



(24) 



where r is the distance between particles, 9 is the angle 
between r and z, 



r = y/Jx — x') 2 + (y — ?/) 2 + (z + z 1 ) 2 

is the distance between particle 1 and image of the par- 
ticle 2 and 9 is the angle between f and vertical line (see 

FigEJ. Then cosO = — ^ z+ f z ^ > (we use coordinate 
system with z\\n and z = at the wall so that particles 
have z < under the wall). 



If both particles are located at the same distance h 
iw t 

T»ff 



below the wall then 6 = 5 . Consider the same orientation 
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Figure 5. Image interpretation of the interaction between two 
particles near the homeotropic wall. This barely supplements 
the exact formulas (11211. 



Formulas (|22p - (|2"4")l can be simply treated, as well, us- 
ing image interpretation. It is clearly seen that the 
interaction between two particles 1 and 2 near the 
homeotropic wall consists of the usual direct interaction 
U\i and image interaction between particle 1 and image 
of the particle 2 ( and vice- versa 2 1; image of the 
particle p has dipole moment —p and image of p' has 
dipole moment —j/j quadrupole moments of the images 
are equal to quadrupole moments of real particles. See 
Figj5|). As each of two image interaction has weight fac- 
tor 1/2 we come to the formulas ([2"2"j) - (f2"4"]l . 



Figure 6. (Color online) Effective power of repulsion between 
two particles at the same height h near the homeotropic wall 
in (12611 vs. dimensionless distance r/h. Red line 1 corresponds 
to the dipole-dipole interaction, first term in (|25|1; black line 
2 - quadrupole-quadrupole interaction, third term in (|25l) . 
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Figure 7. (Color online) Effective power of repulsion between 
two particles at the same height h near the homeotropic wall 
in (I26[) vs. dimensionless distance r/h. Corresponds to the 
sum of all interactions (|25[) for a — 2, j3 = 0.5 and radius of 
particles a\ = a2 = 0.5ft. 



of dipoles p — aa\,p' — era 2 ,, c = —j3a\, c! = —/3a 
total energy of interaction between them: 



\. Then 
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jrwall 

same, horn 

\-nKa\a\ 



8h 2 



(Ah 2 



2)5/2 



6a/3(cii + 02) 



(3r 2 - h 2 )h 



(Ah 2 



2)7/2 



(25) 



+3,5 aia 2 



I20h 2 



560/i 4 



r 5 (4h 2 +r 2 ) 5 / 2 (4h 2 + r 2 ) 7 / 2 (4h 2 + r 2 ) 9 / 2 



r 



This potential is repulsive, elsewhere. Let us present 
it, approximately, in the form of power law dependence 
with some effective power that depends on the distance, 



dipole and quadrupole moment, we show on the FigjT] 
dependence of the effective power on the dimensionless 
distance r/h for the real case a — 2, (3 — 0.5 and particle 
radius a = 0.5h (all lengths should be in terms of the 
height h). 



U, 



wall 

same. horn 



c 



(26) 



_ dlogU 



-XT' — 



where j e ff may be found as j e ff - --g^r 
Fig|5]shows dependence of such effective power on the di- 
mensionless distance r/h for dipole-dipole (a 7^ 0,/3 = 0, 
see red line 1 on the FigJS] ) and quadrupole-quadrupole 
(a = 0,/i ^ 0, see black line 2 on the FigJBJ). We see 
that on small distances r/h <C 1 effective powers are 3 
and 5, respectively. For large distances there is some 
crossover from 3 to 5 for the dipole-dipole interaction, 
with the value near 4 for medium distances. This is sim- 
ilar to electrostatics, where, similarly, dipole interaction 
is weakened from the third to the fifth power near the 
conducting grounded wall. This power 7 e // ~ 4 was ob- 
tained, as well, in computer simulations in [2g|. 

Quadrupole interaction has the same fifth power for 
very large distances r/h S> 1 but in the middle area 
2 < r/h < 5 changes greatly with effective power being 
j e fj ss 5.9 for r/h — 2.5 and 7 e // ~ 3.8 for r/h = 5. 
As real particles with hedgehog configuration have both 

I 



a) 



b) 



c) 



Figure 8. Different configurations of particles near planar 
wall, a) Usual dipole-dipole, b) p-p configuration, c) h-h con- 
figuration according to [7|. 

VI. INTERACTION BETWEEN TWO 
PARTICLES NEAR ONE PLANAR WALL 

In the planar cell axis zllrioo is parallel to the rubbing, 
so that the height of the particles is denoted as a; € (0, 00) 
(see Fig. |8]). Using the Green function (fTS]) for this case, 
we obtain necessary potentials as shown, herein, before: 



jrwall t n I 

U dd , plan = PP {1 _ 3cos 2 e) _3PP 
AttK r 3 y ' f 5 



1 2 • 2 
(p sin tf 



-2 2/ 

r cos t 



(27) 



Tjwall 
dQ 7 plan 

AttK 



, , , COS0 2 , , Ibpcosip 
= (pc — cp ) — — (15cos v — 9) + (pc — cp ) 



•2-2 -2 22 

p sin if — r cos t 



(28) 



jrwall 1 
U QQ,plan _ CC 

AttK ' r 5 1 



90cos 2 6 + 105cos 4 6») + ^ [105 p 2 cos 2 (p - f 2 ){p 2 sin 2 <p 



-2 2/ 

r cos t 



(29) 



r 



Here p = y/(z — z') 2 + (y — y') 2 is the parallel projection 
of the distance r between particles (p lays horizontally 
and makes the angle tp with z direction ). As usually 9 
is the angle between r and z, 



is the distance between particle 1 and the image of the 
particle 2 and 9 is the angle between f and vertical line. 



Then cos9 — 



(x+x 



> and 9 is the angle between r and 



f = \/Jx + x') 2 + (y - y') 2 + (z- z') 2 
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A. Dipole-dipole interaction at the same height for 
the same dipole orientation, FigJHla 

If both particles have the same height h = h' then 
p = r and ip = 9. In this case the dipole-dipole energy of 
interaction (J^TJ) between two parallel dipoles at the same 
height h near the wall with planar anchoring conditions 
(Fig|HJa) can be presented as: 




Figure 9. (Color online) Map of the attraction and repulsion 
zones for dipole-dipole interaction between two particles at 
the same height h near the planar wall. At short distances 
r < h attraction and repulsion cones coincide with unlimited 
nematic case (see Fig[2]a) and make angle 9 — Arccos (l/s/S) 
with z axis. But at great distances r > h repulsion lateral 
cones dramatically collapse from the right and left side into 
dumbbell-shaped region. Sign '-' means attraction,'-!-' means 
repulsion (radial force f r = — 9v dd^pi an ^ f r < o and f r > 
respectively). On the thick blue line f r = 0. This effect is 
absent in usual electrostatics for dipoles near grounded wall. 

Let's analyze dipole-dipole interaction (|3Tj) which cor- 
responds to the usual dipole configuration on the FigJSJa. 
First remarkable behavior is the dipole-dipole isotropic 
attraction on far distances. Really for any finite h we 
have asymptotic form of the dipole-dipole potential for 
great distances r >• h 

2 

Vdd,plan ~ (32) 

FiglHl represents attraction and repulsion zones obtained 
from the potential (|3ip . At short distances r < h attrac- 
tion and repulsion cones coincide with unlimited nematic 




Figure 10. (Color online) Dimensionless energy of dipole- 
dipole interaction vs. dimensionless distance r/h. a) Energy 
for all angles 8 £ (0, n), clear attraction for 8 — 0, it; b) En- 
ergy U± in perpendicular direction 9 = tt/2 between r and 
rubbing direction z. New effect: for r > Vdd = 4.76fo attrac- 
tion appears. 




Figure 11. (Color online) Energy of interaction U± in kT be- 
tween two dipole beads near planar wall in the perpendicular 
direction to the rubbing vs. distance r in pm for different 
heights from the wall: l)/i=2 pm, 2)/i=3 pm, 3)/i=4 pm. 
Here ai = a2 = lpm, a = 2.04, /3 = 0.5. 



case (see Figj2ja) and make angle 9 — Arccos(l / y3) with 
z axis. But at great distances r > h lateral repulsion 
cones dramatically collapse from the right and left side 
into dumbbell-shaped region. On the thick blue line ra- 
dial force is zero f,. = — dVd ^p' an = 0. It crosses perpen- 
dicular direction y/h'm. the point r = = 4. 76 ft.. This 
means that potential Vj^ plan has minimum at the dis- 
tance rdd — 4:.76h, so that for r > rdd there is attraction 
and for r < rdd there is repulsion in the perpendicu- 
lar direction to the rubbing. Actually this minimum is 
unstable and of the saddle type as particles tends to at- 
tract along 9 — 0,7r directions (see Fig fTUj) . This effect 
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of attraction in perpendicular direction along 9 = n/2 
near the planar wall may be found only on far distances 
T > rdd where the interaction is weak. On the FigQj] we 
show the total potential U± in kT units for the beads with 
radius a = lpm which carry dipole and quadrupole mo- 
ments a = 2.04, f3 = 0.5 on the distances between them 
in pm. Beads are located near the wall on three different 
heights h = 2pm, h = 3/im, h = Apm. It is clearly seen 
that effect can be measurable only for h = 2 — 3 pm. Ef- 
fect of dipole attraction on far distances near the planar 
wall was observed as well in (2(| with help of computer 
similations. 

Actually we can say that potential between two dipole 
particles is approximately becomes isotropically attrac- 
tive (|32|) only for distances r > Tdd- This effect is absent 
in usual electrostatics for dipoles near grounded wall. 



B. Dipole-quadrupole interaction at the same 
height for the same dipole orientation, FigJSJa 

If both particles have the same height h = hi then 
p = r and ip = 9 and we can find dipole-quadrupole 
potential from (|28| . Then it can be presented in the 
form: 

jjwall,same 

A dQ rfT 2 = a P( a l ~ a 2 )VdQ,plan(r, h) (33) 

AitKa^a^ 



Map of the attraction and repulsion zones of this po- 
tential is presented on the FigfTJ] (we consider greater 
particle is in the center a\ > a 2 )- At short distances 
r < h attraction and repulsion cones coincide with un- 
limited nematic case (see FigfJJb) and make angle 6 = 
j4rccos(3/-\/l5) with z axis. But at great distances r > h 
lateral attraction and repulsion cones dramatically col- 
lapse from the right and left side into dumbbell-shaped 
region. Sign '-' means attraction,' +' means repulsion (ra- 
dial force f r = - 9Vd ^ lan ; f r < and f r > respec- 
tively). On the thick blue line f r = 0. It crosses axis y/h 
in the point r = rd q = 7.29h so that at great distances 
t > Tdq this potential can be approximately presented in 
the form: 

&COS0 

VdQ,plan (,-JOj 

which is repulsive elsewhere in the upper half-plane z > 
and is attractive elsewhere in the lower half-plane z < 0. 
This effect as well is absent in usual electrostatics for 
dipoles and quadrupoles near grounded wall. 

C. Quadrupole-quadrupole interaction at the same 
height for the same dipole orientation, FigJHJa 

If both particles have the same height h — h' then 
p = r and ip = 9 and we can find quadrupole-quadrupole 
potential from (|2"9l . Then it can be presented in the 
form: 



VdQ,pla 



cos9{l5cos 2 9 - 9) 15cos9r(sin 2 9r 2 - Ah 2 ) 



(4/i 2 + r 2 ) 7 /2 



(34) 



jjwall,same 
U QQ,plan 

AirKa\al 



/3 2 aia 2 VQQ, p iari(r, h) 



(36) 



rubbing 



z/h 




Figure 12. (Color online) Map of the attraction and repulsion 
zones for dipole-quadrupole interaction (|33[) between two par- 
ticles at the same height h near the planar wall with greater 
particle in the center. 



with 



90 



plan =Ax + A 2 cos 2 9 + A z cos A 9 (37) 



9 (6O/1 2 - 15r 2 



r 5 (Ah 2 + r 2 )V 2 ' 



15r 2 



105r 4 - 420r 2 h 2 



{Ah 2 + r 2 )V 2 (4h 2 +r 2 ) 9 / 2 



105 



105r 



(Ah 2 + r 2 f/ 2 



(38) 



Map of the attraction and repulsion zones of this po- 
tential is presented on the FigfT51 At short distances 
r < h attraction and repulsion cones coincide with un- 
limited nematic case (see FigfJJc). But at great dis- 
tances r S> h repulsion lateral cones from the right 
and left side collapse into dumbbell-shaped region while 
upper and low cones expand from 9\ = Arccos(0. 86) 
to 9[ = Arccos(l/V5) . Sign '-' means attraction,'-!-' 
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nibbing 



z/h 



+P 2 a 1 a 2 V Q Q^i an (v 1 h) 

Antiparallel configuration of h-h type (p = —p' 
-aa 2 , see FiglSJc): 




U 



wallah — h 
plan 



= -a 2 Vdd, P ian{r, h)+a/3(ai+a 2 )VdQ,pian(r, h)- 

(42) 



+P 2 a 1 a 2 VQQ, p i an (r, h) 

In these configurations dipole beads repel each other but 
in p-p case the repulsion is weaker than in the h-h case be- 
cause of dipole-quadrupole interaction. In p-p case there 
is dipole-quadrupole attraction and in h-h case there is 
dipole-quadrupole repulsion. 



Figure 13. (Color online) Map of the attraction and repulsion 
zones for quadrupole-quadrupole interaction (|36[1 between two 
particles at the same height h near the planar wall. 



means repulsion. On the thick blue line radial force 
fr = — 9VQ ^ lan — 0. It crosses axis y/h in the point 
r = r qq = 7.00/i so that at great distances r > r qq this 
potential can be approximately presented in the form: 

-6 + 30cos 2 6 



VQQ.plar. 



(39) 



which is repulsive elsewhere in the upper and low cones 
with < 9 = Arccos(l / y/E) . This effect as well is ab- 
sent in usual electrostatics for quadrupoles near grounded 
wall. 



VII. INTERACTION OF THE ONE PARTICLE 
WITH NEMATIC CELL 

In this section we consider two parallel walls with dis- 
tance L between them filled with nematic liquid crystal. 
In other words this is a usual nematic cell with thick- 
ness L. We shall consider both cases of homeotropic and 
planar nematic orientations in the cell. 



A. Interaction of the one particle with 
homeotropic cell 

Consider first one particle in the nematic homeotropic 
cell. Let h be the distance from the particle to the top of 
the cell. We choose z = at the upper plane of the cell, 
z > above and z < below so that particle has z = —h 
and bottom of the cell has z = —L.(sce Fig fT4|) . 



D. General elastic interaction between two 
particles at the same height, Figj8ja,b,c. 

The total energy of interaction between two dipoles 
with parallel orientation at the same height h near the 
wall with planar anchoring condition (FigJHJa) can be pre- 
sented as: 



U, 



wall, same 
plan 



4nKa 2 a 2 



a 2 V dd ,pian(r, h) + a/3(ai - a 2 )V d Q ,pian{? , h) 

(40) 

+P 2 a 1 a 2 VQQ^i an {v, h) 

For untiparallel configurations we should take p' = —p. 
There are two possible cases. Antiparallel configuration 
of p-p type (p = —p' = aa 2 , see FigfSJb): 



IT. 



wallop— p 



plan 



^a 2 V d d,pian('r, /i)-a/3(ai+a 2 )V" d Q :P ; a „(r, h)-\ 

(41) 



L 
' 2 



I 

I I 
1 I I 



-I 



Figure 14. One particle in the homeotropic cell. 5 is the shift 
from the center of the cell. Equilibrium shift S = does not 
depend on the thickness L of the cell (|48|) . See Fig |15l as well. 
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If we consider all images of the unit charge located in 
the point z = —h then we can construct Green function: 



C/iom( X J X ') — 



fffl(x,x') = - £ 



1 



H%l(x,x') (43) 



1 



n— — oo 



|A»(x)-x'| ^ |B„(x)-x'| 



where A n (x) = (a;, y, —z+2nL) and B n (x) = (x, y, z + 
2nL) for x = (a;, y, z). Substituting into (ITU brings 

energy of the one particle within the homeotropic cell 



Tjhom,cell tj 

U self ~ U dd,self 



horn, cell _j_ jjhom,cell 



u. 



dQ,self w QQ,self • 

— 1 :A,rV 



self, ceil _ nKa a ^ 1 



dd,hom 



^ \nL + h[ 



self , cell _ 
dQ,hom ~ 



sign(nL + ft) 



71— — OO 



|ni + ft| 4 



C(3) 
(44) 

(45) 



self .cell 
QQ,hom 



3nK/3 2 a 



2„6 



E 



2„6 



|nZ + /i| £ 



3tt K {3 a 

Z 1 " 



<(5) 
(46) 



is Riemann zeta function, £(3) 



where = £ n=1 
1.202, C(5) = 1.036. 

Here h is the distance from the particle to the upper 
plane of the cell, '-' corresponds to p-configuration and 
'+' corresponds to ^-configuration, see FigEJ We see 
that dipole-quadrupole energy (|45|) is zero when particle 
is located in the middle of the cell h = 4 . It should be 
zero as well from the symmetry reasons as in the middle 
of the cell there is no difference between configurations 
with p < and p > 0. 

Self energy of the single quadrupole par ticle within the 
cell was obtained as well in the paper [271 ] using the coat 
approach (but authors there did not find dipole-dipole 
and dipole-quadrupole interaction with the cell). Au- 
thors received formula similar to (|46[) but with unknown 

multiplier V and without the constant 37rK / 5 g6 C(5). If we 
set T = 2nKc — —2f3TrKa 3 in [27j we receive the same 
formula f4"o| up to a constant. 

From this expressions we can obtain approximate en- 
ergy of the one particle in the homeotropic nematic cell 
near the center of the cell if we consider first terms of 
series with n = 0,-1. Then for one particle which car- 
ries dipole moment a ^ it's energy in the homeotropic 
nematic cell is: 



Trjself,cell A-rrKa 2 ^ 

U dd,hom ~ mn " U 



1 



1 



{L + 25Y (L-28Y 



± (47) 



±247rifa/3a 5 



1 



1 



(L + 28) 4 (L - 2SY 



Then minimum condition 



brings equilib- 



rium shift of the particle from the center (see FigfT5|): 



Sen = ± — 0, 

a 



(48) 



where '+' corresponds to p-configuration, '-' corresponds 
to /i-configuration (see FigJ3]). This means that this 
shift from the center of the cell is positive (up) for p- 
configuration and negative (down) for /i-configuration. 
Surprisingly that 5 eq does not depend on the thickness L 
of the cell! 

U - U center, kT 
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Figure 15. (Color online) One particle with radius a\ = 2/xm 
in the homeotropic cell with thickness L = lOfim. Depen- 



dence of the energy shift U — U z 



■-L/2 



on the shift 8 from the 



center of the cell. Blue line 1 corresponds to the dipole par- 
ticle a = 2, /3 — 0.5. Equilibrium shift 5 — 0.5 is in line with 
(148 [) . Thick dashed black line 2 corresponds to the quadrupole 
particle a — 0, /3 — 0.5 . 

If dipole moment is zero, then energy of the one 
quadrupole particle in the homeotropic nematic cell is 
approximately : 



1 



(49) 



(L-26) 5 , 

where S is the shift of the particle from the center of the 
cell z = (5 eq = in equilibrium). 

B. Interaction of the one particle with planar cell 

In this subsection we consider one colloidal particle 
within the nematic cell with director parallel to the 
planes. Let the rubbing direction be along zllrioo. We 
choose coordinate system for this planar cell with axis x 
looks down < x < L and L is width of the cell (see 
Fig 113 ) so that particle has the height x = h from the 
upper plane. 

Using the same procedure as in the Sec. Ill, B (see for- 
mula (US])) we construct Green function for the planar 
cell from the Green function (|4"3")) for the homeotropic 
cell: 
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Figure 16. One particle in the planar cell. 8 is the shift from 
the center of the cell. Equilibrium shift 8 — 0. See Fig |17l as 
well. 



^pian( X ! X ') 



+ J ff^„(x,x') 



(50) 



11 plan\^i A I 



oo 

E 

n=—oo 



i 



A„(x)-x'| 



OO 

E 



i 



|B„(x)-x'| 



where A„(x) = (— x + 2nL,y,z) and B„(x) = (x + 
2nL, y, z) for x = (x, y, z). 

Substitution into (|TT|) brings energy of the 

one particle within the planar nematic cell U 



j jplan.cell . jrplan,cell 
U dd,self ' U dQ,self 



V, 



plan, cell 
QQ,self : 



u. 



self, cell _ S-n-Ka 2 ^ 
dd,plan j 



self 



2„4 



^ \nL + h\ 3 



71— — OO 



nKa a 
2L 3 



U, 



self .cell 



15nK(3 2 a 



u: 



2„6 



self, cell 
an 



<(3) 
(51) 
(52) 



QQ 7 plan 



1G 



E 



97tK/3 a' 



2„6 



|ni + /i| E 



8i 5 



C(5) 
. (53) 

Self energy of the single quadrupole particle within the 
cell was obtained as well in the paper (27j (but authors 
there did not find dipole-dipole and dipole-quadrupole 
interaction with the cell) . Authors received formula sim- 
ilar to (|53|) but with unknown multiplier T and without 

the constant 97r ^ a<> C(5). If we set T = -2(3ttKcl 3 in 
[27| we receive the same formula (|5"5)) up to a constant. 

Expressions (|5H53I) show that the equilibrium position 
of the particle within the planar cell is the center of the 
cell z = -§. We can obtain similarly approximate energy 
of the one particle in the planar nematic cell near the 
center of the cell if we consider first terms of series with 
n = 0,-1. Then for one particle which carries dipole 
moment it's energy in the planar cell: 



U, 



self , cell 



dd,plan 



6nKa a 



2„4 



1 



(L + 2(5) 3 (L-2Sy 



(54) 



If dipole moment is zero, then energy of the one 
quadrupole particle in the planar nematic cell is approx- 
imately: 



j rself ,cell 
qq,plan 



2„6 



1 



(L + 25) 5 (L-25Y 



(55) 



where S similar is the shift of the particle from the center 
of the cell z = h. 

In the paper [l5[ authors used formula (|47|) obtained 
for homeotropic cell to find dependence of the shift 5 
from the cell thickness L in the planar cell. This is not 
correct though they received some qualitative agreement 
with experimental data. In fact there formula (l54l) for 
the planar cell should be used to fit experimental results 
more accurately. 



U - U center, kT 




Figure 17. (Color online) One particle with radius a\ = 2\im 
in the planar cell with thickness L = 10/im. Dependence of 



the energy shift U — U z 



=L/2 



on the shift 8 from the center 



of the cell. Blue line 1 corresponds to the dipole particle 
a = 2, /3 = 0.5. Thick dashed black line 2 corresponds to the 
quadrupole particle a = 0, ft = 0.5 . 



VIII. INTERACTION BETWEEN TWO 
PARTICLES WITHIN THE HOMEOTROPIC 
NEMATIC CELL 



In this section we shall consider interaction between 
two particles located within the nematic cell with 
homeotropic alignment on both walls. The thickness of 
the cell is supposed to be L (see Fig (pi 



Actually we already have Green function (|43[) for this 
case. But it is much more convenient to use another 
form of this Green function which is commonly used in 
electrodynamics (see f32l|): 
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. oo oo , 

^(x,x') = l-£ £ e-(^') si ^ sin ^7 m (^<)X m (^>) (56) 

n— 1 m— — oo 



Here heights z,z', horizontal projections p <1 p > and 
Im,K m are modified Bessel functions. Then using of (fT2)) 
brings dipole-dipole interaction in the cell : 



U; 



horn 
dd,cell 



l6irKpp' 



(nn) sin— —sin — = — Kq(—t—) 
±j ±j ±j 

n=l 



Dipole-quadrupole interaction is: 



(57) 



jjhom 
u dQ,cell 



-^}J(^) 3 ^ (-^)x (56 



n=l 



pc • sin— —cos — — 
Li ij 



■p c • cos—j—sm- 



We see that if particles are located in the middle of the 
cell z = z 1 = jfr then dipole-quadrupole interaction is zero 
^dQcell = 0- Similar quadrupole-quadrupole interaction 
takes the form: 



jjhom 
U QQ,cell 



IQttKcc' 



L" 1 



n=i 



n7rz n7rz ,nnp 

'COS — —COS — ; — Kq( ) 



4 1M« 

(mr) cos c 



L 



with p being the horizontal projection of the distance 
between particles. FigfTS] demonstrates application of 
this formula ([59| for the repulsion potential between two 
spherical particles (with planar anchoring on the sur- 
face providing quadrupole director configuration a = 
0, (3 7^ 0) with diameter D = 2a = 4Apm in the 
center of homeotropic cell (z = L/2) with thicknesses 
L = h = 6.5pm and L — h — 8pm ( experimental 
data are taken from [l6[ ). The only unknown param- 
eter j3 — 0.28 ( we remind that c = ~f3a 3 ) fits both 
thicknesses very well in the energy scale 1 -j- lOOOfcT. 
Spherical particles in this experiment [l6| have boojum 
director configuration (see FigfTJ:.) and there is no exact 
estimations for j3 for this case. But the value ft = 0.28 
is very reasonable and is in accordance with experimen- 
tal estimation from [6J, where authors found j3 — 0.2 for 
hedgehog configuration FigfTJ a. 

The similar formula (|59[) was obtained in the paper 
28] within the coat approach with unknown multiplier T 
as the parameters of the coat remain unknown and it is 
difficult to estimate them correctly for the case of strong 
anchoring condition Wa/K > 1. If we set T = 2wKc = 
—2(3irKa 3 in [28| we receive the same formula (1ST 



IX. INTERACTION BETWEEN TWO 
PARTICLES WITHIN THE PLANAR NEMATIC 
CELL 

In this section we shall consider interaction between 
two particles located within the nematic cell with planar 



Z 
L0 



I T 

|l I 
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- x 



Figure 18. Two particles in the homeotropic cell with thick- 
ness L. 



h = 1.8 x D 

h = 1.5 x D 




x/D 



Figure 19. (Color online) Experimental data taken from [16j 
- energy of elastic interaction between two spherical particles 
with diameter D — 2a = 4Apm in the homeotropic cell with 
thicknesses h = 6.5pm and h = 8pm. Solid blue and green 
lines are taken from the formula (|59|) . The only unknown 
parameter ft = 0.28 fits both thicknesses very well in the 
energy scale 1 ~ lOOOfcT. 



alignment on both walls. The thickness of the cell is 
supposed to be L ( see Figl20l ). 

Although we already have Green function (|50|) for 
this case it is much more convenient to use another 
form of this Green function which can be obtained from 
([56| with help of coordinate system transformation like 
in Sec III,B. Let's turn coordinate system (CS) of the 
homeotropic cell CS hom (x, y, z) round the y axis on ir/2. 
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in 



Then we will have CS plan (x,y,z) with transition matrix 
A: x = Ax, x' = Ax' so that x — z,y = y, z = —x. Then 
Ghom(x,x') = G> iom (Ax, Ax) = Gpi a „(x,x). Omitting 
sign ~ we may write Green function for planar cell in the 
CS plan with n\\z and x perpendicular to the cell plane 
(x€[0,L]): 



Figure 20. Two particles in the planar cell with thickness L. 



oo oo 



n— 1 m— — oo 



w v ' sin sin I m ( )Km\ ) 

L L L i 



(60) 



where heights of particles x, x', horizontal projections p< = \J y 2 + z 2 , p> = ^/ y' 2 + z' 2 , tgy) = ^,tgcp' — p- and p< 
is less than p> . Then taking derivatives from (|12p brings all interactions in the planar cell. Dipole-dipole interaction: 



U. 



plan 



dd,cell 



1? 



{Fx - F 2 cos 2 if>) 



(61) 



where 



00 / \2 / 

E(nir) nirx nirx 
sin sin 
9 T, T, 



71=1 



K o{—)+ K 2{ — )\ -{nir) cos—^-cos—^—K ( — ), 



OO 

F 2 = y^(n7r)^gm : ^gm ::: ^X2( ::: 7^)- 
Lj lj lj 

71=1 



with p being the horizontal projection of the distance 
between particles and tp is the azimuthal angle between 
p and z. 

Dipole-quadrupole interaction: 

U dQ™eii = ~JjT( pc ' ~ c P') cosi P{Ci + C 2 cos 2 Lp) (62) 
where 

C 1 =L [ F ip - 2 Jl 



c 2 = l(^-f^ p j, 

Quadrupole-quadrupole interaction: 



U, 



plan 



WttKcc' 



QQ.cell jj, 

where 



{D 1 + D 2 cos 2 <p + D 3 cos 4 ip) (63) 



D 2 = L 2 - 
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When both particle 
center of the cell x 

Fx 
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L 
2 

— [KoC-T) ^ 

and F 2 



are 

x 

(mr) 2 



111 

we 



the 
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3cos 2 tp) for p <C L. In the limit of 
big distances p > L we have = 1 — — + o( — ) with 

accuracy 5% already for p = L. So for p > L the radial 
component of the force between particles may be written 
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Figure 21. (Color online) Map of the attraction and repulsion 
zones for dipole-dipole interaction (|61|l between two particles 
at the center of the planar cell with thickness L. Director 
no 1 1 z. At short distances attraction and repulsion cones co- 
incide with unlimited nematic case (see Fig[2]a) and make 
angle 9 = Arccos(l/y/3) with z axis. But at great distances 

2 

r > L attraction cones are deformed to the parabola z — . 

2 

Black thin line is the parabola z — -g-. Sign '-' means at- 
traction,'-)-' means repulsion. On the thick blue line radial 

av plan 

force f r = - T ce " = 0. 

as f„ = = -±<^F 2 '>) ■ (1 - ^ - o»V) 

so that dipole-dipole interaction is attractive (f p < 0) 
for -y> c < <p < ip c , ip c = arccos{Jl - ^) « and 

is repulsive for ip c < (p < 2tt — ip c (if dipoles are parallel 
each other p = p' and vice versa if p = —p' ). In other 
words for p > L dipole-dipole interaction is attractive 

2 

inside parabola z — and is repulsive outside this 
parabola (see Fig |21~|) . 

Similarly in the limit of small distance p <C L be- 
tween particles we come to the result of unlimited ne- 
matic U^ n cell = *jj£(pd - cp')costp(15cos 2 <p - 9) and 

U QQ n ceii = ^^(9 - 90cos 2 p + 105cosV) for p « L. 

Map of the attraction and repulsion zones of the dipole- 
quadrupole potential is presented on the Fig |22l (we con- 
sider greater particle is in the center a% > CI2). At short 
distances r < 0.5L attraction and repulsion cones coin- 
cide with unlimited nematic case (see FigJ5]b) and make 
angle 9 = Arccos(3/ vl5) with z axis. But at great 
distances r > 0.5L upper and low cones shrink to the 
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Figure 22. (Color online) Map of the attraction and repulsion 
zones for dipole-quadrupole interaction (|62[) between two par- 
ticles at the center of the planar cell with thickness L with 
greater particle in the center. Black thin line is the parabola 
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Figure 23. (Color online) Map of the attraction and repulsion 
zones for quadrupole-quadrupole interaction (163[) between two 

particles at the center of the planar cell with thickness L. 

2 2 

Black thin lines are the parabolas z — 22_ anc l y = ir- 
respectively. 

2 

parabola z — 

Map of the attraction and repulsion zones of the 
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quadrupole-quadrupole potential is presented on the 
FigHSl At short distances r < L attraction and repulsion 
cones coincide with unlimited nematic case (see Fig|2]c). 
But at great distances r > L upper and low repulsion 

2 

cones shrink to the parabola z = and lateral repul- 

2 

sion cones shrink to the parabola y — ^jh. 



X. DISCUSSION 

Here we want to discuss the current approach in the 
context of others theoretical approaches used for descrip- 
tion of colloidal particles in NLC. We will emphasize ini- 
tial assumptions of them as well as relationships between 
them. 

1. The current approach ([2(], [3(| and this paper). 
Initial assumptions. 

Let one particle creates director field on far distances 
n M (r) = p^ + 3c^§ with /i = x,y. This means that only 
things that other particles 'see' and 'feel' are these two 
constants p and c through the director field. Therefore as 
initial points of this approch we take these two constants 
plus principle of least free energy. So this approach can 
be called: constants plus principle of least free energy. 
Nothing else. 

Then it enables to introduce some effective functional 
F e f f (jHl) which gives necessary solutions on the extremal 
of this functional. And this potential gives all necessary 
potentials between particles as a consequence (TTTj) . (TT2"1) . 

2. Coat approach [H, [51 H3, HI- This approach 
was suggested by authors in [231 ] . 

Initial assumptions of the coat approach are the follow- 
ing: some volume of the NLC around the particle which 
contains all strong deformations and defects is called the 
coat. Outside the coat deformations are small. On the 
surface of the coat some anchoring coefficient W(a) is ar- 
tificially introduced (1st assumption) which depends on 
the point a of the coat so that surface energy takes the 



And now it is apparent that quantitative estimates 
of these constants a^^, (3 3su , 73^ obtained in [23| are 
wrong though the form of the functional is true (as we 
think) . In the case of axially symmetric particles we have 
ct^n = and the functional takes the form ©. It is nec- 
essary to abandon simple connection between constants 



form: 

F s = <£ do-W{a)(nvf (64) 

J coat 

Symmetry of the W(a) coincides with the symmetry of 
the director field around the particle. Plus possibility 
of the gradient expansion of the director is supposed 
(2nd assumption) adn^ < 1. So the coat approach = 
gradient expansion of the director adn^ -C 1 plus un- 
known anchoring W(o~). In the case of small particles 
without defects the coat coincides with the particle itself 
W(cr) — W sur f ace . But even in this case one assumption 
is necessary: coat approach —gradient expansion of the 
director adn^ <C 1. 

Then it enables to find interaction potentials for any 
form of the coat and to find the connection between the 
symmetry and the type of the interaction potential [23j j . 
As well this approach was used for description of the 
spherical particle in confined NLC - in homeotropic and 
planar nematic cell [271 |28[ . Surprisingly this approach 
gave the same results as in formulas (|2"Tj) . (|4*6")l and ([55)1 
but with unknown multiplier F ! If we set F = 2ttKc = 
-2/3nKa 3 in [27ll28| we receive the same formulas (|21l) . 
(|46l) up to a constant and ([59]) . 

But gradient expansion demand adn^ <C 1 can not be 
satisfied quantitatively on the surface of the particle as 
director is changed on the distances comparable to the 
size a of the particle. Really for dipole case %i(r) = p^ 
we have adn^ ss apfo? ~ a ~ 1 near the particle. There- 
for the coat approach can not give exact quantitative 
results. Nevertheless it is very surprising that it gives 
similar results as (|2"I]). (gBJ) and ([55]). 

Apparently the reason of such resemblance is the fol- 
lowing: in the coat approach another effective functional 
is obtained as intermediate result after gradient expan- 
sion on the coat surface. Actually there surface term 
takes the form up to a constant [23j : 

F s = a^rifj, + (l^dsu^ + 73 S fe M 9 s 9 fc n^ (65) 

with constants , /3 3sL , , 73 S fe M taken as surface integrals 
on the coat. They play the role of charge, dipole and 
quadrupole moments. So that total effective functional 
obtained there has the form: 



(66) 

I 

and the surface and to consider constants as just con- 
stants. As far as this connection is not trivial and can not 
be obtained within the bounds of suggestion adn^ <C 1. 
It takes to make experiment or to make computer simu- 
lation to find these constants. Then approaches 1. and 
2. are in fact equivalent. 



F C oat = K I d 3 xl + a^n^ + hsudsn^ + j 3 sk^d s d k n 
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3. Approach of JH, [Hj. 

In this approach the particle is surrounded by imag- 
inary sphere which contains all defects inside and the 
director is supposed to be fixed firmly on the surface of 
the sphere. Outside the sphere director follows Laplace 
equation. Thus this fixed director on the surface plays 
the role of the source of deformations in the bulk NLC. 
Then Green function on the sphere enables to find di- 
rector field in the whole space and to find interaction 
potentials between particles. So authors do not make 
gradient expansion as in (23|. 

Qualitatively results of this approach conforms with 
two previous approaches. Similarly interaction potential 
between particles was obtained as the multipolc series 
expansion [24[ and even three particles effects have been 
found. This approach was applied as well for the inter- 
action between one dipole particle and the wall. Com- 
parison of (|15p and (TK)1) shows that repulsion of the 
one particle with dipole moment from the planar wall 
is 3/2 times stronger than from the homeotropic wall 
U^:ff a Vu'Z l ff aU = 3/2. Authors of |H as well ob- 
tained this ratio to be 3/2. 

Nevertheless this approach differs quantitatively from 
the previous ones. Approach of [24], [25| predicts dipole- 
dipole force to be three times weaker and quadrupolc- 
quadrupole five times weaker than results of [2C| for 
the same director field on far distances. This discrep- 
ancy has no roots in different definitions for dipole or 
quadrupole moment. If someone consider director field 
in [2~ij near the particle and rewrite it in the form 
n^r) = p^ + 3c^f with \i = x,y and then compare 
elastic interaction between two particles in the same def- 
initions p and c then dipole-dipole force will be 3 times 
weaker and quadrupole-quadrupole force will be 5 times 
weaker than results of (2Cj for unbounded NLC. 

Actually this discrepancy is unclear now. 

All three approaches contain some additional assump- 
tions for the simplification of the problem. No one of 
them solves the problem in the exact formulation de- 
scribed in Sec. II . 

Until now all experimental efforts about colloidal par- 



ticles in nematic liquid crystals were concentrated on the 
study of interaction potentials and structures of the par- 
ticles only [l[-[l6|. But experiments on the precise mea- 
surement of the director field near one particle are absent 
absolutely. Exactly this kind of experiments which can 
measure precisely director field and elastic interaction be- 
tween particles simultaneously can be a testing area for 
validity of the one approach or another. 

XI. CONCLUSION 

To conclude we apply method developed in [30l | for the- 
oretical investigation of colloidal elastic interactions be- 
tween axially symmetric particles in the confined nematic 
liquid crystal (NLC) near one wall and in the nematic cell 
with thickness L. Both cases of homeotropic and planar 
director orientations are considered. Particularly dipole- 
dipole, dipole-quadrupole and quadrupole-quadrupole in- 
teractions of the one particle with the wall and within the 
nematic cell are found as well as correspondent two par- 
ticle elastic interactions. Set of new results has been pre- 
dicted: effective power of repulsion between two dipole 
particles at height h near the homeotropic wall is reduced 
gradually from inverse 3 to 5 with increasing of dimen- 
sionless distance r/h; near the planar wall the effect of 
dipole-dipole isotropic attraction is predicted on the large 
distances r > rdd = 4.76/i; maps of attraction and repul- 
sion zones are crucially changed for all interactions near 
the planar wall and in the planar cell; one dipole particle 
in the homeotropic nematic cell was found to be shifted 
on the distance 8 eq from the center of the cell indepen- 
dent of the thickness L of the cell. The proposed theory 
fits very well experimental data for the confinement effect 
of elastic interaction between spheres in the homeotropic 
cell taken from [16| in the range 1 -j- 1000/cT. 

These surface induced effects show the differences be- 
tween nematostatics and electrostatics which can be 
tested experimentally. All numerical calculations in the 
paper were performed using Mathematica 6 and in all 
series we used summation Xm=i ■ 
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